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ABSTRACT 

Several models for type la-like supernovae events rely on the production of a self-sustained detona- 
tion powered by nuclear reactions. In the absence of hydrogen, the fuel that powers these detonations 
typically consists of either pure helium (He) or a mixture of carbon and oxygen (C/O). Studies that 
systematically determine the conditions required to initiate detonations in C/O material exist, but 
until now no analogous investigation of degenerate He matter has been conducted. We perform 
one-dimensional reactive hydrodynamical simulations at a variety of initial density and temperature 
combinations and find critical length scales for the initiation of He detonations that range between 1 
- 10 10 cm. These sizes are consistently smaller than the corresponding Chapman- Jouguet (CJ) length 
scales by a factor of ~100, providing opportunities for thermonuclear explosions in a wider range of 
low mass white dwarfs (WDs) than previously thought possible. We find that virialized WDs with 
as little mass as 0.24 M can be detonated, and that even less massive WDs can be detonated if a 
sizable fraction of their mass is raised to a higher adiabat. That the initiation length is exceeded by 
the CJ length implies that certain systems may not reach nuclear statistical equilibrium within the 
time it takes a detonation to traverse the object. In support of this hypothesis, we demonstrate that 
incomplete burning will occur in the majority of He WD detonations and that 44 Ti, rather than 56 Ni, 
is the predominant burning product for many of these events. We anticipate that a measure of the 
quantity of 44 Ti and 56 Ni produced in a helium-rich thermonuclear explosion can potentially be used 
to constrain the nature of the progenitor system. 

Subject headings: White dwarfs — Nuclear reactions, nucleosynthesis, abundances — Novae, cata- 
clysmic variables — Supernovae: general 



1. INTRODUCTION 

The study of type la supernovae (SN la) occupies a 
critical position with regard to our understanding of the 
Universe. SN la have already delivered powerful insights 
into several areas of physics, especially cosmology and 
nucleosynthesis. Yet, despite their immense value to as- 
trophysics, how SN la actually explode remains a vex- 
ing problem. It has been generally accepted that these 
violent explosions are the result of the thermonuclear 
disruption of compact carbon and oxygen (C/O) white 
dwarf (WD) stars that have neared, or exceeded, the 
Chandrasekhar mass limit. Beyond that, however, little 
is certain. For instance, the mechanism(s) by which SN 
la progenitors reach their critical states is currently un- 
known, and more fundamentally, we remain ignorant as 
to what these critical conditions are. Complicating the 
issue further, WDs only need to be close to the Chan- 
drasekhar limit to be capable of auto-igniting at their 
centers, and external mechanisms can lead to the ignition 
of even less massive WDs. These mechanisms include 
collisions between WDs in dense stellar environments 
or in multiple stellar systems (Rosswog et al. 2009a; 
Raskin et al. 2009; Hawley et al. 2012; Katz & Dong 
2012), tidal encounters with moderately massive black 
holes (Luminet & Pichon 1989; Rosswog et al. 2008; 
Ramirez-Ruiz & Rosswog 2009; Rosswog et al. 2009b; 
Haas et al. 2012), mergers of WDs with neutron stars 
or stellar mass black holes (Fryer & Woosley 1998; 
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Fryer et al. 1999a,b; Lee & Ramirez-Ruiz 2007; Metzger 
2012; Fernandez & Metzger 2013), or the accretion of 
dense material from a WD companion (Nomoto 1982; 
Woosley et al. 1986; Bildsten et al. 2007; Fink et al. 
2007, 2010; Guillochon et al. 2010; Shen et al. 2010; 
Woosley & Kasen 2011; Dan et al. 2011, 2012; Sim et al. 
2012). SN la initiation and the diversity of possible out- 
comes once ignition occurs are therefore pressing issues 
to the astrophysics community that call for further study. 

Several models have been proposed for the sequence of 
events that lead to a SN la-like events, but all share the 
common feature that a self-sustained detonation powered 
by nuclear reactions propagates through a large fraction 
of the interior of the host object (or objects). In the 
absence of hydrogen, the fuel that powers this detona- 
tion consists of either pure helium (He) or the prod- 
ucts of He burning, typically a mixture of carbon and 
oxygen. This detonation must be triggered by an in- 
crease in density, temperature, or both. The conditions 
for the initiation of the detonation is a topic that has 
been discussed at length in the context of C/O burning 
(e.g. Seitenzahl et al. 2009), but only to a limited extent 
for the ignition of pure He (Khokhlov & Ergma 1986; 
Khokhlov 1989). In this paper, we present a parame- 
ter space study to find the critical temperature gradient 
length scales necessary to successfully initiate a detona- 
tion in He media for a given combination of density and 
temperature. 

A detonation will occur when the nuclear timescale is 
on the order of the local dynamical timescale. However, 
the definition of both of these timescales is ambiguous; 
the effective nuclear timescale depends on a number of 
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Fig. 1. — The p — T Parameter Space Relevant to He Ignition. 
The blue-shaded region encompasses the p — T values for which 
a He detonation might occur, and is defined by the constraint 
0.1 < Td yn / r 3u < 10. The z/-filled region above approximately 
eight billion Kelvin shows where the energy losses due to neutrino 
emission (pair annihilation, plasma, photoneutrino, recombination, 
and bremsstrahlung processes, Itoh et al. 1996) dominate over the 
energy generation due to triple-o; reactions (Khokhlov &; Ergma 
1986; Rosswog et al. 2008). The red-dashed line represents the 
p — T trajectory of virialized n = 3/2 polytropes (Hydrostatic He 
WDs) up to O.45M0 in total mass; the core density and tempera- 
ture for a O.45M0 WD are denoted by the solid red dot. The green- 
dashed lines show the time evolution of double white dwarf (DWD) 
mergers with a O.4M0 accretor and donors of 0.8Mq, 0.9Mq, 
I.OMq, and I.IMq, where the plotted line thickness increases with 
donor mass, which were produced using the simulations presented 
in Dan et al. (2012). The orange-dashed line corresponds to the 
head-on collision of two OAMq He WDs, and were obtained from 
the simulations of Rosswog et al. (2009a). The purple-dashed line 
exhibits the parameters of a O.2M0 WD as it is tidally disrupted by 
a 5OOOM BH at /3 = 5, and is taken from Rosswog et al. (2009b). 
The dashed cyan line shows the trajectory of the accretion disc 
formed by the merger of a 0.3Mq WD and a 1AM@ NS calculated 
in Metzger (2012). Arrows are used throughout to indicate time 
evolution. 

reactions with different rates, and the effective dynamical 
timescale depends on the size of the region being heated. 
Furthermore, the heated region is most likely not uniform 
in either composition or temperature at any particular 
time, and even less so over the timescale in which the 
heating occurs. We can analytically estimate the length 
scale £ that must be heated to produce a detonation by 
calculating the speed at which the detonation will travel 
u and then multiplying by the timescale of the corre- 
sponding nuclear reactions r nuc , £ = ur nuc . Both u and 
T nuc can be calculated by using simple one-dimensional 
models of the detonation in combination with a nuclear 
network (Khokhlov 1989). This estimate is an upper 
limit on the size required, as heating a region larger than 
this scale would result in an overdriven detonation that 
would eventually slow to the steadily propagating so- 
lution. However, while this method can establish the 
length scale of the resulting detonation structure, it is 
unclear if a region smaller than £ may be capable of yield- 
ing a detonation, especially given that the initial condi- 
tions are unlikely to resemble the assumed steady-state 
detonation structure. For instance, the hydrodynami- 
cal simulations of He detonations in WD atmospheres 
of Townsley et al. (2012) show that steadily propagating 
weak detonations can often strengthen towards a steady- 
state solution, though only after the detonation has trav- 



eled some distance. 

To determine the critical temperature gradient length 
scale for a given density and floor temperature, while 
simultaneously resolving the initiation region and sub- 
sequent shock wave front, we perform one-dimensional 
reactive hydrodynamical calculations that begin with a 
hot spot in a He medium and check if they lead to deto- 
nations. Our measurement of these critical length scales 
allow further constraints to be placed on SN la models 
as well as other more exotic thermonuclear explosions in- 
volving the disruption of He WDs. For many systems, 
we find that the spatial scales relevant for the initia- 
tion of these He detonations are unresolved in even the 
most computationally expensive multidimensional sim- 
ulations, with a continuum of critical lengths £ cr i t that 
ranges from 1 cm to 10 10 cm. 

This paper is structured as follows. In Section 2, we 
present the SN la initiation models in question and the 
density and temperature conditions in which they sub- 
sist. In Section 3 we describe the numerical methods 
and initial conditions used for the study and present the 
results. In Section 4 we discuss the results and their 
application to thermonuclear explosion models. 

2. RELEVANT CONDITIONS 

White dwarfs, the end point in the evolution of 
the vast majority of stars in the universe, are ex- 
tremely common, accounting for about 65% of the to- 
tal heavy element abundance (i.e. heavier than He, 
Fukugita & Peebles 2004). Since WDs experience no 
mass loss through single stellar evolution, their ma- 
terial can only be ejected through dynamical interac- 
tions with other objects and/or by the rapid release 
of nuclear binding energy. Although at present the 
reservoir is dominated by C/O WDs, which possess a 
thin layer of He (Lawlor & MacDonald 2006), a non- 
negligible fraction of the mass lies within pure He WDs 
(with masses in the range of O.17-O.45M , ~6% of 
WDs, Brown et al. 2010; Kleinman et al. 2013) and hy- 
brid WDs consisting of small C/O cores and thick He en- 
velopes with total masses of ~O.5M (Iben & Tutukov 
1985; Nelemans et al. 2001b). Pure He WDs cannot be 
the result of the evolution of single stars within the life- 
time of our Galaxy; these WDs can only be produced 
via interacting binary systems (Nelemans et al. 200 lb, a; 
Rappaport et al. 2009) and approximately half of all 
WDs in close binaries are expected to be either He-core 
or hybrid WDs (Ruiter et al. 2011). In addition to sys- 
tems in which both stars are WDs, binary systems in 
which the accretor is either a neutron star or a black 
hole are also possible. As a result, the conditions under 
which He may be ignited are diverse, with a variety of 
accretion rates, initial entropies, and surface gravities. 

In Figure 1 we present the density-temperature pa- 
rameter space relevant to He ignition. If the timescale 
on which the He material can react, its dynamical 
timescale Td yn = (Gp) -1 / 2 s, is much shorter than the 
burning timescale r^ a — 9.0 x 10 -4 Xg exp (4.4/Tg)/?g 2 
s (Khokhlov & Ergma 1986), the material can expand 
rapidly enough to quench burning (T = 10 9 Tg K and 
p = 10 6 p6 g cm -3 ). This is because in a pure He environ- 
ment the combustion rate is limited by the triple-a reac- 
tion. Appreciable burning will therefore only take place if 
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T3a ^ Tdyn . This comparison of timescales can be used as 
a simple estimate for whether or not substantial He burn- 
ing can occur in particular physical systems, but the on- 
set of dynamical burning does not guarantee that a self- 
sustained detonation will form, as the conditions for det- 
onation are non-trivial (Khokhlov 1989). However, we do 
not expect that detonations will occur in systems where 
T3a ^> Tdyn, since burning would likely be quenched by 
expansion. Above temperatures of approximately 8 x 10 9 
K, energy loss by neutrino emission begins to overtake 
the energy production due to triple-a reactions (the gray- 
shaded region in Figure 1). We have therefore restricted 
our simulations to the region in which rs a < 10rd y n 
and in which the relation \E V \ < E^ a holds, where E v 
includes contributions from pair annihilation, plasma, 
photoneutrino, recombination, and bremsstrahlung pro- 
cesses. For comparison, we have included trajectories for 
virial hydrostatic He WDs (n = 3/2 polytropes with total 
mass of up to O.45M ), double WD mergers, collisions 
of WDs, tidal compression of WDs orbiting black holes, 
and the disruption of He WDs orbiting neutron stars in 
Figure 1. These trajectories all fall within the region of 
interest, and thus are suitable candidates for systems in 
which successful He detonations can occur. 

3. CONDITION FOR SUCCESSFUL HELIUM 
DETONATIONS 

3.1. Numerical Methods and Initial Conditions 

To study the detonation of He, we have carried out 
a set of ID cartesian simulations using the adaptive 
mesh refinement code Mezcal. The code has been 
tested extensively against standard test problems, and 
has been used to simulate magnetohydrodynamic jets 
(e.g., De Colle & Raga 2006; De Colle et al. 2008) and 
gamma-ray burst afterglows (De Colle et al. 2012b, a). In 
Mezcal, the equations of hydrodynamics are integrated 
in time using the second-order Runge-Kutta method, 
while the fluxes are calculated using the HLL method, 
with a second-order spatial reconstruction of the primi- 
tive variables at the cell interfaces. A thirteen element 
"a-chain" nuclear network which includes all alpha pro- 
cess elements from 4 He to 56 Ni (Timmes et al. 2000) is 
integrated by an unsplit solver, and includes losses due 
to neutrino cooling. We employ the Helmholtz equation 
of state (Timmes & Swesty 2000), which includes radia- 
tion, ion, and electron-positron degeneracy pressures. 

To properly couple the hydrodynamics with energy 
injected by burning, and removed by neutrinos, the 
timestep is limited such that the energy change is less 
than 0.5% of the internal energy for all cells. This 
is achieved by setting the time step At™ uc = 5 x 
10 -3 At n-1 eo/|ei — eo|, where eo and e\ are the thermal 
energies before and after the update from the nuclear 
routine. Otherwise, the timestep is set by the Courant 
condition, defined as the minimum among all the cells 
of Ax/(v + c s ), where Ax, v and c s are the size of the 
cell, the velocity, and the sound speed respectively. As 
in other hydrodynamical studies, both the energy release 
calculated by a reduced network and the fraction of 56 Ni 
produced can vary by ~30% when compared to a 489 
isotope network (Timmes et al. 2000), or when using a 
coupling constant that is larger than what is used in this 
study (Hawley et al. 2012). Therefore, while we expect 



that this reduced network should capture the general nu- 
cleosynthetic trends, the normalization of these trends 
may differ from what is presented here. 

For our study we initialized linear temperature gradi- 
ents ("hot spots") of a particular size £ in a pure He 
environment of density p, where the maximum tempera- 
ture T max occurs at x = 0. We assume one-dimensional 
cartesian geometry with reflecting boundary conditions 
at x = and outflow boundary conditions at x = x max . 
For consistency, we always chose the temperature floor 
To = 10 _2 T max . We employed the method of bisection 
to determine the minimum size £ cr i t of the hot spot that 
resulted in a self-sustained detonation for a given T max 
and p. 

As a test of the implementation of the nuclear net- 
work module in the Mezcal code, we first reproduced 
the resolution study for C/O detonations as presented 
in section 3.3 of Seitenzahl et al. (2009), and also did 
not find convergence until the minimum grid cell size 
Ax m in < 2 cm. We then performed a resolution study 
for a pure He composition, with T max = 10 9 K, To = 10 T 
K, and p = 10 6 g cm -3 , and found convergence for 
Ax m in < 10 4 cm, which is approximately equal to the 
estimate of the width of the burning region presented in 
Table I of Khokhlov & Ergma (1986). 

We determined which simulations successfully pro- 
duced detonations by identifying shocks that propagated 
at a near-constant velocity, indicating that the were pow- 
ered by the continual release of energy through sustained 
nuclear reactions. In Figure 2 we show the results of two 
different simulations, where the initial conditions differ 
only in the extent £ of their hot spots. In the left column 
of Figure 2, the temperature gradient quickly evolves into 
a fully- formed shock wave, rapidly converting the major- 
ity of the He fuel into Ni (bottom two panels). The right 
column shows a case where £ has been reduced by a fac- 
tor of ~2. In this case, the burning that is initiated is 
unable to outpace the hydrodynamic expansion of the 
fluid, even though triple-a reactions initially inject a sig- 
nificant amount of energy into the system. The boundary 
between successful and unsuccessful detonation is some- 
times signified by the inability to produce a shock. This 
is in contrast to what Seitenzahl et al. (2009) found for 
the initiation of C/O detonations; they found a shock is 
produced even in sub-critical conditions. 

In Figure 3, we show the time evolution of T and Xm 
for a successful detonation. For p > 5 x 10 6 g cm -3 it 
is common for the initial conditions to lead to multiple 
runaways that only produce a stably-propagating deto- 
nation through their mutual interactions. These multiple 
runaways come as the result of the strong density depen- 
dence of the triple-a reaction. The expansion of the hot 
spot introduces upstream perturbations to the density, 
which can potentially experience dynamical burning if 
the temperature in these regions is ~10 9 K. Each run- 
away produces a pair of left- and right-traveling shock 
waves that interact with the hot spot as it expands, and 
collide both with shocks produced from other runaway 
regions, and with the reflecting boundary condition at 
x = 0. While this behavior is somewhat pathological, 
our criteria for the boundary between success and failure 
is still applicable, as we ensure that these multiple inter- 
actions eventually produce a single shock that propagates 
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Fig. 2. — A Successful Detonation Versus A Failed Detonation. Depicted on the left is the time evolution of a successful (critical size) 
He detonation, initialized with T ma x — 0.25 x 10 9 K, p = 10 6 g cm 6 , and £ cr it = 0.11 x 10 cm. The hot-spot quickly develops into a full 
blown shock-driven detonation. The shock couples with the rapid fusion of approximately 90% of the initial He fuel into nickel, producing 
peak temperatures of ~4 x 10 9 K. On the right is displayed a failed detonation in which the size has been reduced to £ = 0.48 x 10 9 cm, 
a reduction of roughly a factor of two from the critical size; the temperature and density are unchanged from the left panel. This initial 
temperature gradient is unable to form a shock. Despite the initially positive energy generation rate, the products of the triple-cn reactions 
(collectively defined as Ash) are nearly negligible. In both panels the initial state is indicated in red, intermediate states in green, and 
latest states in blue. 
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Fig. 3. — Time Evolution of Temperature and Nickel Fraction 
For A Successful Detonation. The above panels display a successful 
detonation initialized with T ma x = 10 9 K, p = 5x 10 6 g cm 6 , and 
R = 6.1 x 10 4 cm. On the left it is shown that the temperature 
reaches a peak of T = 5.06 x 10 9 K for a brief moment during 
the ignition of the shock wave. The detonation-front subsequently 
relaxes to a temperature of approximately T w 4 x 10 9 K. Three 
ignition phases occur before the shock obtains stability; multiple 
ignitions are characteristic of high density (p$ ^> 1) detonations. 
The ignitions produce shocks in both directions; reflections of the 
backward traveling shocks due to the left boundary are clearly 
visible. In the right panel the nickel production is shown. The 
ignition phases of the detonation convert up to nearly 85% of the 
He fuel into nickel, while the stable shock settles to a conversion of 
about 50%. The large abundance of nickel in the hot-spot region 
is not a general result within the parameter space (c.f. Fig. 7). 



The choice to use a linear temperature profile is admit- 
tedly arbitrary, so we determined the critical size for a 
gaussian temperature distribution and compared the re- 
sults to a linear temperature profile with the same T max , 
To and p. Figure 4 shows the time-evolution of a critical 
detonation for these two cases where £ = £ cr i t (for the lin- 
ear case) or R = R cv \t (for the gaussian case). While the 
delay between the start of the simulation and the time 
at which a self-propagating detonation forms is different 
between the two simulations, the resulting detonation 
structure is almost identical. 

3.2. Critical size estimates 

Past efforts in determining minimal length scales for 
the initiation of detonations have been concerned with 
C/O WDs (Arnett & Livne 1994; Niemeyer & Woosley 
1997; Ropke et al. 2007; Seitenzahl et al. 2009). In all of 
the aforementioned works, these critical sizes were deter- 
mined through the use of one-dimensional hydrodynam- 
ical simulations, as multidimensional parameter space 
searches are computationally prohibitive. But while 
these studies have been performed by numerous groups 
for C/O mixtures, they have never been performed for 
pure He. 

In the left panel of Figure 5, we show the minimum size 
scale £ cr it for combinations of p and T max that success- 
fully resulted in a detonation. These scales range from 
~1 cm at ~10 9 g cm -3 to 10 10 cm at 10 5 g cm -3 , with 
a relatively weak dependence on T at each p. Overlayed 
on this panel are the same regions highlighted in Figure 
1, where the red line again shows the virial relation for 



Helium Detonations 



5 




x [10 7 cm] 

Fig. 4. — Linear Versus Gaussian Initial Temperature Gradients. 
The above detonations are initialized with T max = 10 9 K and p = 
10 6 g cm -3 . The top panel displays a linear hot spot with critical 
radius £ cr it = 7.0 x 10° cm, and for comparison, the bottom panel 
displays a gaussian hot spot with critical radius (defined by the 
half width at half maximum of the distribution) R = 4.3 x 10 6 cm. 
It is evident that the only appreciable difference between the two 
is that the detonation with the gaussian profile reaches the ignition 
point earlier than the linear profile. The times displayed are not 
uniformly distributed. 

WDs. Within the blue crescent (where dynamical burn- 
ing occurs) and along the WD relation, the minimum 
scale ranges from 10 8 - 10 10 cm, similar to the range of 
WD sizes. As shown in Figure 5, the critical T necessary 
to produce a detonation at a given size £ is comparable to 
the virial temperature of the WDs themselves, and thus 
these conditions can potentially be achieved through dy- 
namical processes that are driven by the WD's own grav- 
ity (e.g. collisions or the acquisition of material from a 
companion). 

It is apparent from the left panel of Figure 5 that £ cr it 
is primarily a function of density, with only a slight de- 
pendence on temperature as T approaches the v cooling 
limit and where the dynamical and burning timescales 
are comparable. In the right panel of Figure 5 we show 
how the initiation scale varies as a function of density, 
where the temperature of the upstream material To is as- 
sumed to be 10 7 K, and compare it to the the Chapman- 
Jouguet (CJ) length Iqj = Drs a , where D is the CJ 
velocity for that particular density and r^ a is calculated 
using the downstream temperature (Khokhlov 1989). We 
find that the initiation scales at these densities are such 
that Zqj > 100£ cr it, implying that burning may not pro- 
ceed to completion within the time it takes the detona- 
tion to traverse the system. 

3.3. Nucleosynthesis 

In the simple CJ model of a detonation (Chapman 
1899; Jouguet 1905), the burning is assumed to be instan- 
taneous across an infinitely-thin zone behind the shock, 
and the detonation propagates at a single speed, the CJ 
velocity D. In reality, the reactions that take place oc- 
cur in a finite timespan, and thus the region over which 
the reactions complete has a finite size. Zel'dovich, von 
Neumann, and Doring (ZND) (Zel'dovich 1940; Doring 
1943) improved this by modeling a detonation as a one- 
dimensional leading shock wave moving at the detona- 
tion speed trailed by a reaction zone of finite width 



(Fickett & Davis 1979). This enables one to estimate a 
variety of length scales, including the region over which 
the bulk of the energy is injected, the distance over which 
the fuel is exhausted, and the distance over which the 
composition approaches nuclear statistical equilibrium 
(NSE) (Timmes & Niemeyer 2000). 

For systems in which the CJ and ZND scales are much 
smaller than the scale of the He fuel region, the burn- 
ing is expected to yield a composition of material com- 
posed primarily of He and nickel, since the timescale for 
a-capture is short compared to the triple-a timescale 
(Khokhlov 1989). If the burning is quenched before the 
detonation has traveled a few CJ lengths, the burning 
cannot proceed to completion, as is the case for systems 
in which the initiation scale is comparable to the lo- 
cal fuel scale. Alternatively, the material can expand, 
thereby quenching the burning before alpha particles 
have accumulated on each seed nucleon to produce 56 Ni 
(Townsley et al. 2012). 

Our results are consistent with this picture. In Fig- 
ure 6 we calculate the fraction of 56 Ni in the region im- 
mediately downstream to the detonation within our sim- 
ulations shortly after the detonation reaches a steadily- 
propagating state (left panel). This is compared to the 
expected fraction of 56 Ni if the detonation is allowed to 
propagate a distance long enough for the downstream 
material to reach NSE (right panel). We find that the 
two models produce radically different yields of 56 Ni, in- 
dicating that the post-detonation state that is produced 
shortly after a detonation first forms is not representa- 
tive of the state it would have reached if it was allowed 
to expand until NSE was achieved. 

As described in Section 3.2, the initiation size £ cr it is 
comparable to the size of the WD itself for p < 10 6 g 
cm -3 . Incidentally, this is the same region in which the 
fraction of 56 Ni declines to the point that it is no longer 
the predominant species within the ash (Figure 7, right 
panel). For densities 10 5 < p < 10 6 g cm -3 , and for tem- 
peratures typical of adiabatically compressed WD mate- 
rial (red dashed line), the predominant species that is 
produced is 44 Ti (Figure 7, left panel). This progression 
to heavier species in the ash as the density increases is 
similar to the result seen when He is burned at constant 
density and temperature (Hansen 1971). 

The stable decay product of 44 Ti is 44 Ca, and the pro- 
duction of this particular isotope has been associated 
with type II SNe. However, the amount of 44 Ti pro- 
duced per type II SN may not be sufficient to explain 
the amount required to produce the ratio of 44 Ca/ 40 Ca 
observed in the Sun (Timmes et al. 1996). Given the fact 
that greater than half of the initial He can be burned to 
44 Ti, He burning can produce vastly larger quantities of 
44 Ti than type II SNe, which on their own may not be 
able to explain observed abundance (Woosley & Weaver 
1994; The et al. 2006). If the rate of WD collisions, dis- 
ruptions, and SN la events involving a He WD are com- 
parable to the type II rate, He burning on low-mass WDs 
may be the primary source of 44 Ti in the universe. 

4. DISCUSSION 

We have presented one-dimensional results of the de- 
termination of critical (smallest) spatial scales for the 
gradient initiation of detonations in He matter using the 
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Fig. 5. — Critical Sizes For He Detonation. The left-hand panel shows the minimum size scale ^ cr it for a linear temperature gradient as a 
function of density p and peak temperature T ma x- The T3 a < ^dyn an d \E V \ < Es a shaded zones, along with the virial WD line, have been 
carried over from Fig. 1. <!; cr it ranges from approximately 1 cm to 10 10 cm, and is weakly dependent on temperature. In the right-hand 
panel, we present a comparison of the critical size £ cr it against the width of the CJ burning zone as a function of density. The dashed black 
line represents the Iqj, which is equal to the product of the burning timescale and the detonation velocity (Khokhlov 1989), and the solid 
blue line represents £ cr it assuming an initial temperature To = 10 7 K. 




Fig. 6. — Simulated versus NSE Nickel Abundances. The left panel displays the average fraction of nickel produced by a critical detonation 
within the zone delimited by 3£ cr it < x < 4£ cr i t when the shock front reaches x = 4£ cr i t . These measured values range from zero to ~0.7. 
In the right panel we have calculated the NSE values of nickel abundance under the corresponding conditions. 
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Fig. 7. — Titanium and Nickel Yields In the Hot Spot Zone. The above panels display the average abundances of titanium (left) and 
nickel (right) within the hot spot zone (0 < x < £ C rit) as a function of p and T max - These measurements were taken when the shock front 
crossed x = 4£ cr i t . 
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Fig. 8. — Critical Masses. The minimal peak temperature re- 
quired to detonate a given mass M = p£crit/^© °f density p is 
displayed. Additionally, we have included the trajectories of in- 
ternal mass as a function of density for n = 3/2 polytropes with 
masses 0.15M©, O.25M , O.35M , and O.45M . 



reactive hydrodynamics code Mezcal. In particular, we 
have shown that the spatial scale £ cr i t required for the 
initiation of a detonation differs significantly from the 
CJ size, specifically, it is approximately two orders of 
magnitude smaller for a wide range of densities. That 
the initiation length £ cr it is greatly exceeded by the CJ 
length implies that certain systems may not reach nu- 
clear statistical equilibrium within the time it takes a 
detonation to traverse the He material. The existence of 
a small region in which a detonation can be initiated does 
not guarantee that a sufficient length of fuel is available 
to ensure NSE is eventually achieved. Therefore incom- 
plete burning that terminates at an element lighter than 
56 Ni will occur in many physical systems in which the 
size of the fuel reservoir is comparable to £ cr it- 

In astrophysical systems for which the state of the fluid 
is initially determined by a balance of pressure and grav- 
itational forces, a compressive event will typically oc- 
cur on the local dynamical timescale, which naturally 
results in the material being heated to the local virial 
temperature. Whether this action leads to a detona- 
tion depends on whether or not the region that is be- 
ing compressed has a scale that is comparable to £ cr i t . 
In systems where there is a prolonged period in which 
He experiences repeated compressions, e.g. mass trans- 
ferring or merging binary systems, the criteria for det- 
onation will be satisfied as soon as the length scale of 
the region being compressed is comparable to the ini- 
tiation scale. These repeated compressions are seen in 
stably-accreting, unstably-accreting, and merging sys- 
tems, and are driven by convective motions and/or fluid 
instabilities (Fryer et al. 1999a; Shen & Bildsten 2009; 
Guillochon et al. 2010; Fernandez & Metzger 2013). 

Within an atmosphere in which the mass of the at- 
mosphere is small compared to the mass of the host ob- 
ject, the density falls off exponentially, so material will 
only have approximately constant density over the scale 
height H. As £ cr i t is a strong function of density (Fig- 
ure 5), the detonation is not expected to be sustained 
as it propagates upwards through such a steep density 
gradient, resulting in a "blowout" (Townsley et al. 2012) 
that quenches burning after propagating in the vertical 
direction a distance ~H. However, the detonation can 
continue to propagate in the horizontal direction, as the 
bulk conditions in the base of the atmosphere of a nearly- 



spherical body do not vary much with the position, re- 
sulting in a continual blowout of partially-burnt material 
as the detonation propagates around the host object. 

We find that for densities typical of WDs, the pre- 
dominant element that is produced from burning is not 
necessarily 56 Ni, and is more commonly an intermediate 
a-chain element such as 44 Ti. Since the detonation is 
quenched in such systems after propagating a distance 
comparable to £ cr i t , the results of our simulations can be 
directly used to predict the nucleosynthetic yields result- 
ing from blowouts at a given initial density p (Figures 6 
and 7). For a candidate SN la resulting from a double- 
detonation model, a measure of the amount of 44 Ti pro- 
duced by the event, in combination with a measure of its 
peak brightness, which can be used to estimate the mass 
of 56 Ni produced (Ruiter et al. 2013) (and thus the mass 
of the C/O core that exploded, Sim et al. 2010), may be 
able to uniquely constrain the nature of the progenitor 
system. 

In some circumstances, such as collisions and tidal dis- 
ruptions, the material's temperature is increased in a sin- 
gle dynamical episode to super- virial temperatures, as 
the source of gravity in these systems is external and 
can be arbitrarily large. In these cases, the scale of the 
system need not be comparable to the initiation scale. 
In Figure 8, we transform the critical sizes presented in 
Figure 5 into the mass M contained within the critical 
region, and compare this to the integrated mass corre- 
sponding to a density p within WDs of different masses. 
We find that the temperature increase is sufficient to pro- 
duce detonations in fully- virialized He WDs with masses 
as low as O.24M , although detonations in lower mass 
WDs are possible if the material is raised significantly 
above the virial temperature, as can be the case in ei- 
ther a tidal disruption or in a collision with a more mas- 
sive (possibly relativistic) compact object. This idea is 
given some credence by the results of hydrodynamical 
simulations illustrating that O.2M He WDs do not deto- 
nate after a head-on, equal-mass collision (Rosswog et al. 
2009a), but do successfully detonate when colliding with 
a more massive WD or after being tidally compressed 
and shock- heated by an IMBH (Rosswog et al. 2009a,b). 
However, because £ cr i t is exceeded by the CJ length by 
a factor of ~100, the initiation scales are rarely much 
smaller than the size of the physical systems themselves, 
and thus incomplete burning is observed even in sys- 
tems which are driven to super-virial temperatures (e.g. 
Rosswog et al. 2008). 

While this work presents the first systematic study 
concerning the length of the temperature gradient neces- 
sary to initiate a detonation in a pure He environment, 
the initial conditions we use are a simplification of the 
physical problem. Therefore, caution should be exer- 
cised in the application of these critical radii, which may 
deviate somewhat from what we have presented when 
the simplifying assumptions we have made are relaxed. 
We do find that the critical radii do not depend sen- 
sitively on the functional form of the temperature pro- 
file, as illustrated in Figure 4, but in reality the den- 
sity, temperature, and/or composition are unlikely to be 
homogenous and smooth over the initiation region, as 
they are assumed to be here. Though we have restricted 
our study to one-dimensional simulations and optimisti- 
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cally employed a 13-species a-chain network, we have 
demonstrated that the variety of possible outcomes of 
He detonation is larger than anticipated, and that these 
outcomes have a diversity of predictive applications to 
upcoming transient searches. 
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oration. We thank F. Timmes for making the equa- 
tion of state and reaction network used in this work 
publicly available. We acknowledge support from the 
David and Lucile Packard Foundation, NSF grant: AST- 
0847563, and the NASA Earth and Space Science Fellow- 
ship (J.G.). 
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